On the existence of the Boltzmann-Grad 
limit for a system of hard smooth spheres^ 

M. Tessarotto a - b and P. Nicolini a ' fe 
a Department of Mathematics and Informatics, University of Trieste, 
Italy, b Consortium of Magneto- fluid- dynamics, University of Trieste, Italy 

'do'- 

O ' (Dated: June 30, 2008) 

o : 

Abstract 

S ■ 

Despite the progress achieved by kinetic theory, its rigorous theoretical foundations still remain 
unsolved to date. This concerns in particular the search of possible exact kinetic equations and, 

m ; 

specifically, the conjecture proposed by Grad (Grad, 1972) and developed in a seminal work by 
Lanford (Lanford, 1974) that kinetic equations - such as the Boltzmann equation for a gas of 
■ classical hard spheres - might result exact in an appropriate asymptotic limit, usually denoted 

as Boltzmann-Grad limit. The Lanford conjecture has actually had a profound influence on the 
scientific community, giving rise to a whole line of original research in kinetic theory and mathe- 

7— I ■ 

' matical physics. Nevertheless, several aspects of the theory remain to be addressed and clarified. 

t> : 

In fact, its validity has been proven for the Boltzmann equation only at most in a weak sense, i.e., 

OO ■ 

if the Boltzmann-Grad limit is defined according to the weak * convergence. While it is doubtful 

\o ■ 

whether the result applies for arbitrary times and for general situations (and in particular more 

OO : 

generally for classical systems of particles interacting via binary forces), it remains completely 



unsolved the issue whether the conjecture might be valid also in a stronger sense (strong Lanford 
conjecture). This paper will point out a physical model providing a counter-example to the strong 
Lanford conjecture, representing a straightforward generalization of the classical model based on 
a gas of hard-smooth spheres. In particular we claim that that the one-particle limit function, 
defined in the sense of the strong Boltzmann-Grad limit, does not generally satisfy the BBGKY 
(or Boltzmann) hierarchy. The result is important for the theoretical foundations of kinetic theory. 

PACS numbers: 47.10.ad,05.20.Dd 
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I. INTRODUCTION: BASIC MOTIVATIONS 



Classical statistical mechanics, and in particular kinetic theory represents, is a sense, one 
of the unsolved problems of classical mechanics. In fact, although the microscopic statistical 
description (MSD) of classical dynamical systems formed by A-body systems is well known, 
a complete knowledge of their solutions is generally not achievable. From the mathematical 
viewpoint it provides an example of axiomatic approach following from first principles and as 
such it must be considered as an 'ab initio' formulation. Two equivalent treatments of MSD 
are known, which are based respectively on the introduction of a phase-space distribution 
function (PSDF) either on the A-body phase-space T N or, respectively, on the 1-particle 
phase-space r\. In the T N -approach the PSDF is the so-called microscopic PSDF f N . It 
follows that J'n obeys the Liouville equation, whose characteristics are simply the phase- 
space trajectories of the same dynamical system, to be identified with a classical A-body 
system |6|. This equation is equivalent to a hierarchy of equations (the so-called BBGKY 
hierarchy) for a suitable set of s-particles distribution functions (/i^), obtained letting 
s — 1, .., N — 1, which are uniquely related to the corresponding PSDF. On the other hand, 
in the Ti— approach the PSDF (the Klimontovich probability density k^ N \ defined in the 
Ti— space) evolves in time by means of the Klimontovich equation [2|. Also for this equation 
the characteristics are just the phase-space trajectories of the A— body system, this time 
- however - projected on the Ti— space. Therefore, in both cases it is actually necessary 
to determine the phase-space trajectories of all the particle. Hence, for classical systems 
characterized by a large number of particles (A ^> 1), the computational complexity (of 
this problem) is expected to prevent, in general, any direct calculation of the time-evolution 
either of the A-body or any of the the s-body distributions. This has justified the constant 
efforts placed so-far for the search of 'reduced' statistical descriptions, of which kinetic 
theory (KT) is just an example. This is intended in order to achieve efficient statistical 
descriptions especially suitable for complex dynamical systems, including both gases and 
plasmas. Precisely, the primary goal of KT is the search of statistical descriptions, either 
exact or in some sense approximate, whereby the whole dynamical system is associated only 
to the one-particle kinetic distribution function (/i) defined on the one-particle phase-space 
r 1; without requiring the knowledge of the dynamics of the whole dynamical system. As a 
consequence in KT-descriptions the evolution equation of the kinetic distribution function, 
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to be denoted as kinetic equation, is necessarily assumed to depend functionally, in some 
suitable sense, only on the same distribution function and the one-particle dynamics. In 
particular, one of the most successful developments of KT is doubtless related to the so- 
called 'ab initio' approaches. These are to be intended (in contrast to heuristic or model 
equations) as the KT's which are obtained deductively - by suitable approximation schemes 
and assumptions - from the corresponding exact MSD. In traditional approaches usually KT 



is obtained adopting the Tn -approach to MSD 



Q, EL [ioi 



However, also the Klimontovich 



method (based on the 1^ —approach) can be used [2|, since it is completely equivalent to 



that based on the Tn -approach |24j]. In all cases KT's have the goal of determining the 
evolution of suitable fluid fluid fields, associated to prescribed fluids, which are expressed 
as velocity moments of the kinetic distribution function fx and satisfy an appropriate set 
of fluid equations, generally not closed, which follow from the relevant kinetic equation. 
'Ab initio' kinetic theories are - however - usually asymptotic in character. Namely, kinetic 
equations are typically satisfied only in an approximate (and asymptotic) sense and in a 
finite time interval, under suitable assumptions. 

A. Asymptotic kinetic theories 

A well-known asymptotic kinetic equation of this type is provided by the Boltzmann ki- 
netic equation for a classical gas formed by N smooth rigid spheres of diameter d (Grad,1958 
[l|), which is obtained from the exact equation of the BBGKY hierarchy for the one-particle 
kinetic distribution, i.e., 

F^^t) f[ N) =d 2 (N-l)Cjf\ (1) 

where F\ and C\ are respectively the free-streaming operator F 1 (r 1 ,v 1 ,t) =t| + vi --^ 
and a suitable collision operator If], 10]. For definiteness, in the remainder we adopt a 
dimensionless notation whereby all relevant functions (in particular, the Newtonian particle 
state xi = (ri, vi), the time t , the particle diameter d and the volume of the configuration 
space V) are considered non-dimensional. The transition from the 1-particle equation (pQ) 
can be obtained by adopting a suitable asymptotic approximation and suitable assumptions 
on the joint probability densities^, [J. These require, in particular, the introduction of 
the so-called rarefied gas ordering (RG ordering ), to be meant both in a global and local 
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sense, for the relevant physical parameters. More precisely, by imposing that e = 1/N 
is an infinitesimal, the particle diameter d, the volume V of the configuration space (fl) 
and the particle mass m; the related global orderings are requiring to satisfy the orderings 
(Grad,1958 J) 



d ~ o(s 1/2 ), 

V ~ o(e°), (2) 

7?(r,t) = 47m(r,t)d 3 /3V ~ o(e 1/2 ). (3) 

The last ordering, in particular, prevents the number density n(r,t) from becoming so large 
that volume fraction rj(r,t) can be locally finite, i.e., of order o(e°). In fact, it is well-known 
that if there results locally rj(r,t) ~ o(e°) particle correlations (in particular two-particle 
correlations) may become non- negligible also on the large scale In fact, these 

correlations, which are not generally expected to decay rapidly in time [4|, can be also long 



range in character |25|]. Instead, in validity of the RG ordering defined above, uniformly in 
phase-space and at least in a finite time interval / = [t Q , ti] , with At = t\ — t a such that 
At ~ o(e°), the following conditions are assumed to be satisfied: 

• Assumption #1 - in T s x I ol , the approximate (i.e., asymptotic) joint probability 
densities f s (e) (for any s G N with s N) are smooth and bounded ordinary 
functions defined in T s x I ol , where T s is the s-particle phase-space; 

• Assumption #2 - the asymptotic factorization condition (AFC) 

f s (e, x 1; ..x s , t) = J] Me, x 4 , t) [1 + 0(t - t )o{e«)] (4) 

i=l,s 

is satisfied identically for any s 6 N such that s/N ~ o(e). Here f\ (e,Xj,t) (for i = 
1, s) is the one-particle probability density which satisfies the asymptotic Boltzmann 
equation 

F 1 (r 1 ,v 1 ,t)/ 1 (e,)=rf 2 iVC 1 / 2 (£,), (5) 
and Q(t — t ) is the Heaviside theta function which vanishes for t = t Q ; 
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If the RG ordering and the previous assumptions hold locally (i.e., in the whole phase- 
space Ti and at least in a finite time interval I \ = jtn,U]), the Boltzmann equation ([5]) 



is expected to be locally valid in the same domain [17|, ll9j, |20| at least in an asymptotic 
sense. This means, introducing an arbitrary monotonic decreasing sequence of infinitesimal 
parameters {e} = {ei > 0, i G N} , that the sequence {fi(e, x^t)} defined in terms of them 
is expected to converge in a weak (asymptotic) sense for e — > 0. In other words the whole 
domain Ti x 7 ol [existence domain of fi(e, xi,i)] : 

• the asymptotic solution f\(e, x 1; t) differs, by an error infinitesimal of order o(e ai ) with 
respect to the exact solution /^ (xi,t), being ct\ is an appropriate strictly positive 
real number. As a consequence, the error Aj\ = j\(e)— f[ N \ while remaining non 
zero, can be taken arbitrarily small; 

• in Ti x I i the Boltzmann kinetic equation differs from the exact one-particle BBGKY 
equation at most by terms of order o(e a2 ), where ot<i is an appropriate real number 
< ct2 < 1 in general different from o.\. 

Even if the rigorous proof of the global validity of the Boltzmann equation for arbitrary 
initial and boundary conditions has yet to be reached, its success in providing extremely 
accurate predictions for the dynamics of rarefied gases and plasmas is well known (see for 



example, Cercignani, 1969 



6j; Frieman, 1974 



B. Boltzmann- Grad limit and the Lanford conjecture 

Basic issues remain to be clarified regarding the rigorous theoretical foundations of KT. 
One such problem - and the one we want to address in this Note - refers in particular to the 
search of possible exact kinetic equations and, specifically, the conjecture suggested originally 
w Grad (Grad, 1972 [3j) and investigated by Lanford in a seminal paper (Lanford, 1974 
7|; see also Frieman, 1974 [8]), that kinetic equations - such as the Boltzmann equation 
for a gas of classical hard spheres - might result exact in an appropriate asymptotic limit, 
denoted as Boltzmann- Grad (B-G) limit. In other words, according to this conjecture, there 
should exist a suitable operator L* (denoted as B-G limit operator) such that the limit 
functions f s = L*fs N ^ of the sequences j/i^ j , to be defined in terms appropriate discrete 
sets {Ni G N} , should result exact solutions of the equation of the Boltzmann hierarchy. 



The B-G limit is customarily intended as the limiting "regime" where the total number of 
particles N goes to infinity, while the configuration-space volume V remains constant, the 
particle diameter d goes to zero in such a way that Nd 2 approaches a finite non-zero constant 
and the average mass density Nm/V = M/V remains finite (Grad, 1972 [3]; Lanford, 1974 
Frieman, 1974 8]), i.e., there results: 

1 

N' 
Nd 2 



:,d,m — » 0, 



V 

M = 



-> fci, 
mN 



(6) 



fa 



where ki (i = 1, 2) are prescribed non-vanishing finite constants. In the case of plasmas 
further analogous requirements must be placed on the total electric charge and current 
carried by each particle species Q, . In addition, the proper definition must be made for 
the limit operator L*. In fact, in order that the sequences j/i^ j converge in some sense it 
is necessary to determine their time evolution. According to Lanford and previous authors 
this can be achieved by constructing and explicit solution of the corresponding equation 



of the BBGKY hierarchy, to be represented explicit 



distribution /. 
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13, 



y 



18 



by a time-series expansion for each 
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20, 



21 



22] As a consequence, 



it was found that L* can be defined in the sense of weak* convergence for the sequence 
j/i^ j [3] ■ The proof of the weak convergence of j/i^ j in this sense was first reached in 
the seminal work of Lanford (Lanford, 1974 |7|) who was able to prove also the local validity 
of the Boltzmann equation in a finite time interval I \ = [t , t±] of amplitude At = t\ = t a 
smaller than 1/5 and under the assumption of factorization at the initial time for the joint- 
particle distribution fa. Even if this result does not suffice to justify possible meaningful 
physical applications, the conjecture has actually had a profound influence on the scientific 
community, giving rise to a whole line of original research in kinetic theory and mathematical 
physics. The work was later extended by other authors to include 2D and 3D and global 



validity for the Boltzmann equation. 



by oti 
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Q, 0, Q, 
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221] However, the validity of the Boltzmann equation for general situations remains dubious. 

A key issue, however, is related to the possible validity of the Lanford conjecture in a 
stronger sense, not just for the Boltzmann equation but also for the BBGKY hierarchy itself, 
as following by suitable definition of the B-G limit operator L* acting on the joint probability 
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densities /i^. In fact, let us assume that L* is defined in the sense of uniform convergence 
in phase-space of the sequences to the strong limit functions f s = L* fg N \ In such 

a case the conjecture can be advanced that the strong limit functions f s belong to same 
functional class of j/i^ j (strong Lanford conjecture). In particular, this means that when 
applying the operator L* term by term to the equation of the BBGKY hierarchy for f[ 
[Eq.©] 

L*F 1 ff ) =L*{d 2 (N-l)C 1 fW}, (7) 

the limit function / x = L* f[ N ^ should satisfy the corresponding equation of the Boltzmann 
hierarchy 

F i h = hCif 2 . (8) 

For the validity of this limit equation it follows that L* should commute with the streaming 
operator F 1; in the sense that it should result identically 

[^FjjfW)^, (9) 

being [L*, FJ= L*F i — F i L*. In the sequel we intend to point out, however, that the 
Lanford conjecture is not generally valid in this sense, i.e., that the limit functions defined 
in the strong B-G limit do not actually belong to the same functional class of the sequences 
j/i^l and in particular to the solutions of the Boltzmann hierarchy. Nevertheless, weak 
convergence in the sense indicated above may still be warranted. In order to prove the 
point, in this paper we intend to propose a counter-example based on the introduction of a 
modified three-dimensional hard-sphere problem. 

II. A COUNTER-EXAMPLE: A MODIFIED 3D HARD-SPHERE SYSTEM 

To prove this point, we consider here a system (Sn) of N partially-impenetrable hard- 
smooth spherical surfaces ('spheres'). The other key element of the proof is the adoption of 



the Klimontovich approach. As indicated elsewhere 2J, |26j this permits to construct an ex- 



act explicit integral representation for the s-particle distribution functions without recurring 



to cumbersome time-series representations 24J. For defmiteness, the system Sn is defined 
by requiring that all particles are alike with diameter d and mass m and are included in 
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a bounded and connected 3D configuration space Q of IR 3 of volume V(Q). The particles 
of Spf can be classified respectively as external and internal, according to the sub-domains 
of the SAT-configuration space to which they belong, denoted respectively as external and 
internal {Q ext and Vt int = Q — VL ext ). It is assumed that the two sub-domains are mutually 
inaccessible, i.e., particles cannot move from Q ext to Vl int or vice versa. As a consequence 
the numbers of internal and external particles (defined by the occupation numbers Ni nt and 
iVea^with N = Ni n t + N ext ) are by assumption constant. External particles are those whose 
inter-particle distances (distance between the centers of the same spheres) is larger than 
(or equal) to d. Two particles are called mutually internal if their inter-particle distance 
is smaller than (or equal) to d. Internal particles are, therefore, those such that there ex- 
ists at least another particle of Sn with which they are mutually internal. It is required 
furthermore that: 1) the occupation numbers Ni nt and N ext are both non-zero; 2) particles 
and the boundary of Q are mutually impenetrable; 3) external particles are impenetrable 
when they collide with another (external or internal) particle; 4) two arbitrary mutually 
internal particles are, by definition, mutually impenetrable (since the intersection between 
their boundaries is always non-empty). Particles can undergo interactions either with the 
boundary (unary interactions) or among themselves (binary interactions), all assumed elas- 
tic. For binary interactions, we distinguish between external and internal collisions. In 
particular external collisions occur when two particles - initially with an inter-particle dis- 
tance large than d - touch each other. Instead, internal collisions are defined only among 
mutually internal particles. The MSD for the system Sn, in analogy to the customary 
hard-smooth sphere problem [l|, Q], can be achieved in principle in an elementary way by 
distinguishing between external and internal subsets of phase-space, either r^r or Ti (see 
discussion above). In particular, for example, Tffl and T™ 1 are respectively the subsets of 
in which particles belong respectively to the external and internal sub-domains of the 
configuration space. In the case of the 1^— phase-space formulation, this implies that the 
Liouville equation must be satisfied identically by the PSDF /jv in both subspaces (T^ 
and rj^*). This leaves, nonetheless, a large freedom in the choice of the initial-boundary 
conditions as well as the functional class of . In particular, due to the arbitrariness of 
/at, it is always possible to invoke the assumptions: Assumption a) the PSDF results 
continuous in the whole set and in particular on the boundary between external and 
internal particles (drffi = ST™*); Assumption (3) fx is a smooth and bounded ordinary 
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real function. The corresponding I\— phase-space formulation for Sn, is obtained by con- 
structing the corresponding Klimontovich probability density. In T\ for external particles it 
reads 

k {N - t} (y,t) = -^- £ Sfr-xtWQfat), (10) 

ext i=l,N ext 

where y = (r, v) is an arbitrary state vector of the one-particle phase space T\. Here Xj(£) = 
Xi(x ,t ,t) [for i = 1,N], denote the phase-space trajectories of the particles of Sn with 
initial conditions Xj(£ ) = x Q .j [for i = l,N]. These trajectories are assumed to be defined 
uniquely in the set Tjy x /, / being a suitable bounded time interval, and Qj(r,t) is the 
function G,(r,£) = 1 — J2j=i n j^i®(.d—\r — Yj(t)\), while Q(x) is the Heaviside step function 
( 1 if x > 

Q(x) = \ As a consequence it follows that in the subset of phase-space Ti for 

{ if x < 0. 

external particles (Tl xt ) the one-particle distribution function when expressed in terms of 
the initial microscopic PSPD reads f[ N) (y,t) = J r dx f( N \x ,t )k( N ^(y,t). Invoking the 
Liouville equation for this can be prove to imply: 

f^(y,t) = ff n M-n N) M, (ii) 

where % N \y,t) = J Fn dxf^t)S(y - x x (t)) and jW(y,f) = (N - 
dxf( N \x.,t)8(y — Xj_(t))G(d — |r — r 2 (t)|). Then the following result can be reached 




Theorem - Non-existence of the strong B-G limit for Sn 

Let us assume that there is at least a finite time interval I \ = [t ,ti] C M such that the 
probability densities fg N ^ (s = 1,2) and their strong limit functions f s = L*fg N ^ ( s = 1,2 
) f[ is strictly positive in Fx x J ol , so that there results uniformly in F± x J ol for Sn '■ 1) 
L*l} N) (y,t)) = 0;2) the strong limit function /i(y,t) reads h{y,t) = L*f^ N \y,t);3) fi(y,t) 
satisfies identically the homogeneous equation 

F 1 f 1 (y,t) = 0. (12) 

As consequence, we conclude that in the B-G limit the strong limit function /i(y, t) does not 
generally satisfy equation (Q and hence neither the corresponding equation of the Boltzmann 
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hierarchy. Hence, at least in the case of the hard-sphere system here considered, the strong 
Lanford conjecture for the BBGKY hierarchy fails. 



III. CONCLUSIONS 



In this paper the issue of the validity of the Lanford conjecture in the sense of the 
strong B-G limit has been investigated. An example case has been formulated based on the 
analysis of a system of partially impenetrable smooth-hard spheres. We have shown that 
if the one-particle limit function fi(y,t) is intended in the sense of the strong B-G limit it 
does not generally belong to the functional class {/i} of the solutions of the one-particle 
limit equation. In other words, in such a case the limit function is neither a solution of 
the corresponding equation of the Boltzmann hierarchy nor - as a main consequence - of the 
Boltzmann equation. This result raises obviously the interesting question whether similar 
conclusions can be reached for the customary smooth-hard sphere system jl[ [f| or to more 



general systems of interacting particles. This problem, to 
the approach here developed, will be discussed elsewhere 



gether with a detailed analysis of 



26]. 
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